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SUMMARY 


A simple, direct procedure is developed for converting frequency- domain 
aerodynamics into indicial aerodynamics. The data required for aerodynamic 
forces in the frequency domain may he obtained frcwn any available (linear) 
theory. The method retains flexibility for the analyst and is based upon the 
particular character of the frequency-domain results. An evaluation of the 
method is made for incompressible, subsonic, and transonic two-dimensional flows. 


I!ITRODUCTIC»l 

For many years, unsteady aerodynamic theories and applications have focused 
primarily on the frequency domain since the aerodynamic calculation is simpli- 
fied if the motion of an airfoil or lifting surface is restricted to be single 
harmonic (refs. 1 and 2). However, for applications to aeroelastic systems with 
feedback control and for aeroelastic systems with structural nonlinear ities, it 
is of considerable value to represent the aerodynamic forces ir the time domain. 

For an aerodynamic theory which is linear in the motion oJ the aeroelastic 
system, there is a fundamental correspondence between the frequency and time 
domains through a Fourier transform pair (refs. 1 to 3) . Such a linear theory 
may still include important physical effects such as shock wave motions in the 
transonic regime which are sometimes, incorrectly, thought of as being invari- 
ably nonlinear. In principle, if the aerodynamic forces are known at a suffi- 
cient number of frequencies, a numerical inversion to the tine-domain represen- 
tation is always possible. Such an inversion is rarely practical, however, 
because the aerodynamic forces are only known at a relatively small number of 
frequencies. Instead, a parameter identification approach is usually employed 
whereby time histories of aerodynamic forces are assumed tc be sums of exponen- 
tials. The time constants and coefficients are chosen to give a best fit to 
the fieguency representation of the aerodynamic forces which has been obtained 
numerically by some appropriate aerodynamic solution procedure. 

Such representations in the time domain date from the early work of Jones 
(ref. 4) and extend to the recent interesting work of Roger (refs. 5 and 6) 
and Vepa (ref. 7). Roger employs a particularly straightforward procedure for 
determining his representations while Vepa uses Fade approximants and a least- 
squares method. Abel (ref. 8) and Dunn (ref. 9) have subsequently improved upon 
these methods. It is the purpose here to develop a simple, systematic procedure 
for time-domain representations which retains flexibility for the analyLt and is 
based upon the particular character of the frequency- domain results. An evalua- 
tion of the method is made for incompressible, subsonic, and transonic two- 
dimensional flows. No difficulty is anticipated in using the method for three- 
dimensional flcsws where results are available for the frequency representation 
of the aerodynamics. Finally, although not emphasized here, the general proce- 
dures jiay be useful when it is desired to convert from time-domain representa- 
tions i:o frequency-domain representations. Such applications might arise when 



finite-difference aerodynamic calculal ions lead directly to time-domain results 
(ref. 10). Marc H. Williams, of Princeton University, piovided the frequency- 
domain data used in tlw compressible-flow examples. 


SYMBOLS 

a^ coefficients of exponential time representation 

b airfoil half-chord 

bj exponents of exponential time representation 

C Theodor sen function 

Cl lift coefficient 

lift coefficient due to heaving 

lift coefficient due to pitching 

Cf^ moment coefficient 

Cm' moment coefficient due to heavina 

n 

C^^ moment coefficient due to pitching 

D denominator in polynomial representation of Theodorsen function 

F real part of Theodorsen function 

G imaginary part of Theodorsen function 

h heaving displacement 

I total number of terms in sum 

1 (-1)^/^; also, index for summation 

k reduced frequency, u\b/XJ 

L lift 

M Mach number 

N numerator in polynomial representation of Theodorsen function 

t time 

U free-stream velocity 

a angle of attack; also, angle of pitch 
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fluid density 


T dimensional t'me, Ut^'b 

Waqner function 

transient aerodynamic functions 
' f requency 

Superscript: 

PT piston tiieory 

Suoscr i pts. 

I imaginary part 

R real pai t 

max maximian 

A bau over a symbol denotes Fourier transform; a dot over a symbol denotes 
derivative with respect to time. 


BASIC APPROACH 

For definiteness, csonsicier s<ime aerodynamic generalized foroe, say 
due to sane step change in a motion variable, say H/U. Thus, 


h/ U * 1 

(T ■> 0) 

(U) 

h , U *• 0 

(T < 0) 

Hb) 


Assume may he represented by 




1 




b. T 

a .e ^ 

1 


" 0) (2a) 


Cl - 0 . 


(T - 0) (2b) 


where the are yet to he determined but it is anticipated that bj < 0 

Taking the Fourier transform of equations (1) and (2), 
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(3) 


Ci, ^ ika^ 

h'U _ {'bi ♦ ik) 

i-1 


where a bar above a quantlcy dektotes Fourier tranafora and k ia the transfora 
variable. Taking the real and iaaginary parts of equation (3), 



I/O 


I 

V 


i-1 


a^k2 


bf t k^ 


(4a) 


iTi »i 


(4b) 


At this point, there are two important questions: 

(1) Is a representation such as equation (3) or equations (4) capable of 
matching the known behavior to arbitrary accuracy by increasing the number of 
terms retained in the series? This question is answered in the affirmative by 
numerical examples end, in the special case of inocmpressible flow, the ana- 
lytical results of Desaarais (ref. 11). 


(2) How can a^,b^ oe determined conveniently, simply, and unambiguously? 
Vepa has suggested a (modified) least-squares procedure for determining aj 
and b£ Here a simpler procedure is used. The are determined by the 

extreaa of then the a^ are determined by a least-squares fit to the 

frequency-domain data for only, subject to the two constraints that 


the real part is identically satisfied at k ■ 0 and The resultant 



is then predicted at intermediate 


k values. 


Moreover it is assimied 


that the bj, which are the poles of the aerodynamic transfer function (see 
eq. (3)), are independent of the particular generalized force and motion and 
are inherent characteristics of the dynamics of the fluid. Vepa anticipated 
this assiaption would be useful but did not pursue it. The procedure is shown 
in this paper to give good results. Dunn (ref. 9) has adopted this assiaption 
partially by using the same b^ for each distinct type of motion, e g. , hea\<e 
and pitch. 
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rdr oaapleteness, Boger's procedure is also briefly described here 
(refs. S and (} . A uxlaw value of reduced frequency w is selected which 
is an upper liaiit on the frequency range of interest. Vest the bj are chosen 

as 



(i - 0, 1, 2, 3, . . I) 


The a£ are then deterained by a least-squares prooedure using both real and 
iaaginary parts of the aerodynanic transfer fimctions (aatrix elements). 

Another characteristic of Roger's procedure, thou^ not absolutely ettsential, 
is that the prooedure is applied to the aerodynasic influence aatrix relating 
pressure to downwash, rather than to the matrix relating generalized farces to 
generalized ccordir^atc This automatically imures all motions and resultant 
aerodynamic forces are treated on a common basis. Finally, the limits k 0 
and k • itre not enforced as oonstraints in Roger's method. Abel (ref. 8) 
has modified Roger's method to enforce the oonstraints at k « 0. 

From both a theoretical and practical point of view, it is better to select 
only the imaginary part of the aerodynamic transfer function to construct the 
representation and to allow the real part to be predicted. From a practical 
point of view, this approach provides an internal check against: (1) numerical 

errors in the frequency-domain data and (2) deviations of aerodynamic data from 
linearity in the motion if they are taken from experinent or finite-difference 
calculations. From a theoretical point of view, tte real and imaginary parts 
are those associated with a single time-dependent function. Thus, oonstructing 
a valid representation of the imaginary part is sufficient to insure a valid 
representation of the real part. In principle, of course, an alternate approach 
is to construct the representation using the real part of the aerodynamic trans- 
fer function and predict the imaginary part. It will become clear in the fol- 
lowing examples, howevv;r, why this alternate approach is not the preferred 
choice. 


H«OOM>RESSraLE FLOW 

For linear, potential, incompressible, two-dimensional flow, the fluid 
unsteadiness is characterized in the frequency domain fay a single function, 
the Theodersm function C(k). It is related to the Wagner function 4> which 
is the lift due to a step diange in downwash at the airfoil three-quarter chord 
(usually said tc be a step change in angle of attack) through a Fourier trans- 
form pair (refs. 1 and 2), i.e., 


♦ <T) 


1 

2ff 


r il!ll 

‘ ik 


dk 


(5) 
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( 6 ) 


c(lf) - ik dT 


Folloving the basic approach, assiaie that ^ may be represented by 
I 

♦(T) • y (T > 0) (7a) 

1-1 

♦(T) -0 (T < 0) (7b) 


Using equations (() and (7), the corresponding representation of Theodorsen's 
function is 


C(k) 


2 (-bi ♦ ik) 
i»l 


( 8 ) 


or, in teras of its real and iaaqinary ooaponents C - P -«■ iC, 



(9a) 


G 



aibik 

2 2 
bf * k 


(9b) 


The frequency-doaain results for F and G are well known (refs. 1 and 2) 

and are shown as dashed lines in figure 1. The qu<*stion is how to deteraine ai 
and b{. First consider the ai. As k 0, F 1; and as k F 1/2. 

These liaits arc well known for any aerodyuaaic theory, since k 0 is the 
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steady-flow Halt and k * ^ is given in general ithough not for M « 0) by 
the piston theory (refs. ! and 2). 

Requiring equation (9a) to satisfy these liaits aeans b-, « G and 
a) ■ F(k»0) « 1 so that 


*1 ^2 ■*■••• * F(k»«») » 1/2 


( 10 ) 


Consider now equation (9b) and, for siinlicity, let I « 2 so that only b 2 
reaains to be determined. There are various ways this might be done, e.g.: 

(1) Collocation, i.e., require b 2 to be sudi that equation (9a) or (9b) 

is exactly satisfied at some intermediate k 

0 < k < • 

(2) Least squares, i.e., require equation (9a) and/or (9b) be satisfied 

in a least-squares sense 

Strictly, this second approach leads in general to a nonlinear equation for 
the bj^. Vepa (ref. 7) avoids this difficulty by rewriting equation (S) as 
a ratio of two polynomials: 


C(k) 


M(ik) 

D(ik) 


( 11 ) 


and then multiplying equation (11) through by D, i.e., 

D(ik)C(k) - M(ik) (12) 


before applying the least-squares procedure to determine the coefficients In the 
polynominals of N and D. This approach does lead to linear algebraic equa- 
tions for these coefficients. While the procedure put forward below is no more 
rigorous fundamentally than that of Vepa, it avoids two objections which might 
be raised about Vepa's procedure. First, the b^ retain their individual 
identity and are not lost in complicated expressions for the polynomial coef- 
ficients of N and D} second, equation (8) , or actually equation (9b) , is 
satisfied in a least-squares sense (to determine a^) and not the modified 
equation (12). Satisfying equation (12) in a least-squares sense gives undue 
weighting to high k values. 

The procedure suggested here for determining b 2 is siaqple. For I > 2, 
since b^ ■ 0, 
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( 13 ) 


C 


a2b2k 

“5 5 

b2 ♦ k 


1 


How G ■ 0 at k * 0,a» (which eq. (11) already satisfies) and has an extreatfli 
at k - 0.2. (See fig. 1.) Equation (13) has an extreMtai, using eleaentary 
calculus, at 


k - ib2 (14) 


Bence, select b 2 • -0.2; the Minus si^i gives the correct si^i of G and 
also provides a stable aerodynaaic systea. The corresponding approxiaants to 
F and G are shown as solid lines in figure 1 along with their exact counter- 
parts. The agreaient is reasonable, though certainly iaperfect. 

The approxiaant can be ii^oved by increasing I. The question then 
becoaes how to deteraine the other b^. If there were several extreaa for G 
(and they were well separated as is typical for the iaaginary parts of aero- 
dynaaic generalized forces when Multiple peaks occur), then a value of b£ 
would be selected to be equal to the -k value at each p^ak. In tlie present 
exaaple, however, since no other extrou are present, additional b£ are 
siaply added on either side of -0.2 to iaprove the approxiaation. In figure 2, 
results are shown for b; « 0, b2 ’* -0.1, b 3 • -0.2, and b^ * -0.4 and in 
figure 3 for b; * 0, b 2 ■ -0.05, b 3 " -0.2, and b 4 ■ 0.6. Note that the bj 
on either side of an extreaua are chosen here by inspection and iteration. This 
siiq>lifies the procedure but presuaably incurs soae loss of accuracy cxmpared to 
deteraining the as part of a least-squares solution procedure. The cor- 

responding aj^ were deterained, after the b^ were selected, by a least- 
squares fit to equation (9b) for a selected nuaber of k values subject to the 
constraints of equation (1 0) . The procedure is standard using Lagrange Multi- 
pliers to invoke the constraints, and the details are oaitted. Up to 27 values 
of k were used, although 15 gave essentially the saae results and as few as 5 
gave reasonable results. 

The representations of figures 2 and 3 arc auch iaproved over those of 
figure 1, with those of figure 3 being soaewhat better than tliose of figure 2. 
They could be improved further tsy increasing I. Bowever, this seeas unneces- 
sary; instead, oaaparlsons with alternative representations are considered. 

The well-knotm representation of Jones (refs. I, 2, and 4) is shown in fig- 
ure 4. This corresponds to 1-3 with b^ -0, b 2 • -0.0455, b 3 » -0.3 in 
the present aodel, but the a^ were not deterained b/ Jones using a least- 
squares procedure. Using Jones' values of the b£ and a least-squares proce- 
dure to deteraine the aj^ changes the details of the representation, although 
it is not noticeably iaproved. (See fig. 5.) Thus, it is concluded that Jones' 
approxiaant is less accurate than the present I > 2 approxiaant; the basic 
reason for thia is the nuaber and choice of poles. 
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In figure 6 the are chosen frna Desaarais* continued fraction repre> 
sentation (ref. 11) of Theodorsen'r function for I • 4 (to be disc; '’ur- 
ther in following sections). The results are coaparable to those .eic-tn ea^ . er 
for 1*4, but are not better as far as one can judge. Finally/ ! ^Iiould 
noted that Vepa (ref. 7) using his 3ote elaborate proceduie has also obtained 
excellent representations of F and <>, 

The quantity b^ is ossiaaed to be leai, although froa a mathematical point 
of view/ complex ocnj^'gate pairs are petmissible. Also, all poles are assiased 
to ' simple ones, e.g. / no double poles. Comfdex Md double poles were inves- 
tigered niaaerically/ but their inclusion gar^' no noticeable improvement. This 
result is consistent with Desmarais* continued fraction representation which 
shows that only simple poles exist ed.Ok.^ the negative/ imaginary k-plane axis to 
any order of approximation. Also/ see the discussion ot Edwards in reference 3. 


Continued Fraction Representation 

Let us now turn to a txrief review of the very interesting results of 
Desmarais (ref. 11} for the support they lend to the admittedly heuristic 
procedure described ^bove. Desmarais has established the following continued 
fraction representation of Theodorsen's function: 


C(k) ■ 1 ♦ 


(-1/2) 


1 + 


1 


i4k + 1 


1 + 


i4k 3 


1 + 


i4k + 5 


1 + 


This infinite fraction may be truncated to obtain approximations of various 
orders/ and Desmarais has developed convenient recursion formulas for these. 

At any order of approximation/ C(k) is represented by a ratio of polynomials. 
All of the poles are a).ong the negative imaginary k-axis (corresponding to 
negative/ real b^)/ i.e./ the branch cut of Theodorsen's function. The poles 
become infinitely dense as the order of the approximation is increased. The 
continued fraction representation converges everywhere in the complex k-plane 
except along the branch cut. Thus/ Theodorsen's function has no poles, except 
possibly along the branch cut. 
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the practical ieanoe of the above results is that/ although in fact 

there ace fw poles of C(k)/ one say expect to obtain an approxiaatior. of 
any desired accuracy by representing C(k) as a rational function whoso poles 
are all along the negative, iaaginary k>axis. 

For additional discussion of the continued fraction aodel, see Oesaarais 
(ref. 11). even though there is no known counterpart for coapressible flow. It 
is possible a siailar situation exists. For subsonic flow, at saall k the 
aerodynaaic forces behave very auch as for inooapressible flow, while for large 
k the aerodynaaic forces will aayaptotically approach those of "piston theory.” 
Finally, it is worth aentioning that, for a given niaaber of poles the least- 
squares procedure aay be used to obtain a better aatch with the true C(k) than 
the continued fraction representation (ref. 11). This is not to say that the 
latter is. in general, inferior to the foraer. Indeed, just the opposite is 
true, as w^ll be clear to the reader. It is siaply to say that for the purpose 
of providing sr$ accurate representation of C(k) by a ratio of polynoaials, a 
better niaierical fit can be obtained using the least-squares procedure than that 
given by my specific order trunca' ion of the continued fraction representation. 
It is the fact that the latter aay be used to generate a systasatic and converg- 
ing ressesentation of any order, as well as the theoretical support it gives for 
the least squares aethod, which underscores its fundaaental iaportance. 


Wagner Ftsiction 

Once the a^/b^ are known, one has a representation of the Wagner func- 
tion froa equations (7). Results were obtained for the I « 2 approxiaant (of 
fig. 1) and the better of the two 1-4 approxiaants (of fig. 3). For the 
1-4 representations whose results are shown in figure 3, the corresponding 
Wagner function representation is indistinguishable froa published results 
(refs. 1 and 2). By contrast, the 1-2 result is soaewhat different and it 
is shown along with the exact result in figure 7. Figure 7 and its oo^>ress- 
ible counterparts are the ultimate results of the present procedure. The excel- 
lent agreeaent between the 1-4 approxiaant and the exact result is very 
satisfying. 


APPLICATION TO COMPRESS IB FLOW 

Now let us turn to the effects of ooapressibility. Two exaaples will 
be considered. They are a two-diaentional flat plate at M • 0.7 and an 
NACA 64A00( airfoil at M - 0.84. The niaierical data are, respectively, froa 
clajsical aerodynaaic theory, which assiaes an infinitesiaal perturbation about 
a uniform aean flow, and the transonic aerodynaaic theory of Williaas (refs. 12 
and 13), which considers an infinitesiaal dynaaic perturbation about a non- 
uniform mean flow including a shock wave. The shock motion due to the airfoil 
motion is also taken into account in Williaas* tlieory and is consistently 
treated as infinitesiaal. 
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Flat Plate at M » 0.7 


For the flat plate at M « 0.7, the lift and monent (about midchord) 
due to heaving and pitching (about nidchord) are considered. Hence, four 
aerodynaaic transfer functions are computed: 

These are shown in figures 8, 9, 10, and 11 along with an eight- term (pole) 
representation using the i»:esent -sethod. The were chosen from an examina- 
tion of the imaginary parts of general, it 

is found that the extrena of the imaginary parts of these functions occur at 
the same k vtlues. However, they are more distinct for some functions than 
for others. Pot e,:asiple, compare figure 8(b) to j 

ure 9(b). The latter actually offers a better definition of the extrema than 
the former. The b£ are as follows: 


b] X 0 
b2 > 0.03 

b 3 X -0.1 (where imaginary part has one extremum) 
b4 X -0.3 

bj X -0.8 (where imaginary part has a near extremian) 

bg X -1.2 

bj X -1.75 (another extremun) 
bo * -3.5 


The corresponding are given in table I. The results shown are good repre- 

sentations and no others were studied; however, moderate changes in the b^ 
values and even a reduction in their total number would probably still lead to 
satisfactory results. These were, in fact, suggested 'y/ the results for 

the second example, vdiich chronologically were obtained first. It should be 
noted here that the results for at high k are even better than indi- 

cated as the dominant piston theory term has been subtracted out. See subse- 
quent discussion following the next exanqile. 


Using the above results, the aerodynamic indicial functions were computed. 
These are shown in figures 12 and 13. The definitions of the various aero- 
dynamic terms are: 


h/U ratio of heaving velocity to free-stream velocity 

Cl angle of pitch 


n 


Cl s L/pD^b 

Cu = H/po2b2 

I 

i^ere !• and M ere dimensional lift and moment, respectivelyf about midc^ord 
and b is the dimenEioned hcilf-chord. As with the Inconprescible oounttsrpart/ 
the Wagner function? figures 12 and 13 are the ultimate results of the piesent 
method. They are the time histories the experimentalist ^lould measure and the 
inputs to the ?etoelastic an's equations of motion. At short times (correspcmd- 
ing to high f regencies) the results are e^qiected to be less accurate/ evm\ 
though at T ^ 0 the results ate exact because the enforcement of the piston 
theory constraint in the frequency domain as k ^ 

Thr* previously published values (ref. 1) of indicia!, aerodynamic functions 
for pitch and mcsient about the quarter-chord are available in a somewhat dif- 
ferent form* For heaving, the nondimensonal lift and moment are defined ty 




I. 

Pd 2 h 
2n -^(2b)- 


<!>( 


M 


2n -^(2b)2hU 


For pitchirg, the published results are for a step change in pitching velocity^ 
but zero pitch angle - a mathematically well defined but physicsQ.ly turtificial 
motion. No compzurisons were made for this case since comparisons of th' results 
for heaving which are shown in figure 14 were so encouraging. Note that not all 
of the differences between the present results and those previously published 
should necessarily be attributed to inaccuracies in the present approach. See 
Ashley's discussion on pages 347 to 350 of reference 1, as well eic the original 
papers cited in reference 1. Also see Edwards (ref. 14). 

Vepa (ref. 7) has made a similar comparison to published results using his 
procedure for M = 0.5. Simileu: agreenent (and differences) were noted. 


raCA 64A006 Airfoil at M = 0.84 

The example of an NACA 64A006 airfoil at K = 0.84 motivated the present 
work and> in fact, was completed first. Hence, it is considered in somewhat 
more detail, including a study of the effects of number of terms retained in 
the e:qponenti 2 d. time-history representation and, also, the number of k values 
used in the least-squares determination of the coefficients of the ejqKtnentieils. 
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A. R. Seebass,^ of t:he University of Arizona, has also noted the effectiveness 
of such representations for transonic aerodynamics. 

In figures 15 to 18 the transfer functions are shown along with represen- 
tations obtained by the present method. Consider figure 15 first. Results are 
shown for eight and four term representations; and, for the latter, 17 and 28 
k values are used for the least-square' approximation. These results give an 
indication of the sensitivity of the metnod to changes in these parameters. The 
b| are the same as those used in the M « 0.7 flat-plate example, and the aj. 
are given in table II. 

It is clear from <x>mparing the results of figures 15 to 18 to each other 
and to the earlier results for M » 0 and 0.7 that 

(1 ) Hore terms are required at the higher M to obtain a good 
representation . 

(2) Hore terms are required for pitching than heaving motion. 

(3) More terms are required for moment than for lift. 

These conclusions are intuitive but, nevertheless, important. 

Consider in particular the results of figure 18 for moment due to pitching. 
There is a substantial degradation of the representation at high k values as 
the number of terms retained in the representation is reduced from eight to 
four. The effect is exaggerated, however, because in figure 18 (as in fig. 11 
for M » 0.7) the pi s ton- theory contribution (refs. 1 and 2) which is dominant 
at high frequencies has been subtracted out. The piston- theory contribution is 
given by 


_PT U ^'' ab^ 

M \3M ' U 


For a step change in ot, this gives a delta function at i - 0 which is sup- 
pressed in the present presentation of the results. For simple harmonic Kiotion, 
this gives 


PT 4 . 

C = - — (01 «= 

M M 


which clearly dominates for high k over the residual shown in figure 18. 


^ In private communication with the author. 
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Finally, consider the indicial functions which are displayed in £ig> 
ures 1 9 to 22. They are shown for both a four-term and eight-term approxi- 
mation using 28 k values. Comparing these results to those for M «= 0.7 
(cf. figs. 12 13), it is seen the indicial lift behaves in a rather similar 

fashion. However, the indicial moments are different zmd generally smaller in 
magnitude. This fact may explain the relativoly greater difficulty of obtaining 
a good representation of the moments at M = 0.84 compared to M = 0.7. There 
are, of course, no previously pulilished results to which those of figures 19 
to 22 may be compared. 

Using his procedure fc>r supersonic and transonic flow, Vepa (ref. 7} main- 
tains "for higher order approximants the poles behaved in an erratic manner, 
often moving into the right half of the s|s ik] plane." Although Vepa sug- 
gested a possible way of overcoming these problems, the present procedure by 
its nature avoids the difficulty. Soger (ref. 6} has also noted this behavior 
in his work luid chooses the b|^ to be negative to avoid the problem. Sunn 
(ref. 9} allows both aj^ and bj to be determined optimally in a least-squares 
series but does invoke the constraints that tlie bf be negative. 


CONCLUDING SEHASKS 

A simple, direct procedure is suggested for converting freguenc:y-domain 
aerodyna:nics into indicial aerodynamics. The time-domain presentation is in 
the form of a sum of exponentials. Exan^les for classical inconqpressible and 
subsonic, conpressible flow suggest that known results czin be reproduced accu- 
rately. New results are presented for transonic flow based upon Williams' 
frequency-domain theory (AIAA J. , vol. 18, no. 6, June I960). All examples 
studied are two-dimensional; however, no difficulty is expected in treating 
three-dimensional flows where the appropriate frequency-domain aerodynOTlc 
representations are available. 

Not unexpectedly, it is shown that more terms are required in the repre- 
sentation for (1 ) higher transonic Mach numbers (though presumably for suffi- 
ciently high Mach numbers this trend reverses), (2) pitching compared to heaving 
motion, and (3) moment ocnnpared to lift. 


Langley Research Center 

National Aeronautics and Space Administration 
H^unpton, VA 236o5 
August 14, 1 980 
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(b) Imaginary part. 

Figure 6.- Theodor ssn function for I ■ 4, ■ 0, b 2 ■ -0.0 
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14.- Indicial lift and nooent for heaving of a flat plate at N 
(See table I for values of a; and b:.) 
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Figure 18.- Monent due to pitching for £ui NACA 64A00 

values of and bj^. 
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Figure 19.- Indicial lift due to step change in h/U. 
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Figure 20.- Indicial lift due to step change in 



■9 



Figure 21.- Indicial roonent due to step change in h/U. 





